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Multifragmentation scenarios, as predicted by antisymmetrized molecular dynamics (AMD) or 
momentum-dependent stochastic mean-field (BGBD) calculations are compared. While in the 
BGBD case fragment emission is clearly linked to the spinodal decomposition mechanism, i.e. to 
mean-field instabilities, in AMD many-body correlations have a stronger impact on the fragmenta- 
tion dynamics and clusters start to appear at earlier times. As a consequence, fragments are formed 
on shorter time scales in AMD, on about equal footing of light particle pre-equilibrium emission. 
Conversely, in BGBD pre-equilibrium and fragment emissions happen on different time scales and 
are related to different mechanisms. 



PACS numbers: 24.10.Pa; 25.70.Pq 



I. INTRODUCTION 

During the last decade, multifragmentation, i.e. the 
break-up of excited nuclear systems into many pieces, 
has been extensively investigated in heavy ion collisions 
(HIC) around Fermi energies, both from the experimen- 
tal and theoretical points of view [ISHSHIlHEIi]. 
In particular, the study of the mechanism responsible for 
fragment production has driven much attention, espe- 
cially in relation to the possibilityto observe a liquid-gas 
phase transition in nuclei [EE Li El Li- Many efforts have 
been devoted to the characterization of the properties of 
the fragmenting source, such as temperature and density, 
to determine its location inside the nuclear matter phase 
diagram. 

Due to compression and/or thermal effects, the com- 
posite systems formed in HIC may reach low density val- 
ues, attaining the co-existence zone of the nuclear matter 
phase diagram. For instance, an excited system that ex- 
pands under the conditions of thermal equilibrium could 
perform a phase transition staying close to the liquid 
branch of the co-existence line 0, Lul- However, due 
to the Coulomb instabilities, the limiting temperature, 
that a nucleus can sustain as a compact c onfig uration, 
may be lower than the critical temperature (Tj. In this 
situation, the system is brought inside the co-existence 
zone of the nuclear matter phase diagram and undergoes 
a spontaneous phase separation, breaking up into several 
fragments 0,0. 

It is generally believed that in central HIC at Fermi 
energies, the composite matter can be compressed up to 
twice the normal density value (as revealed for instance 
by the emission of energetic particles through hard two- 
body scattering [lj), and then the system expands and 
breaks up into many pieces 0, IE • O ne °f the most 
challenging and still open questions is the understanding 
of the fragmentation mechanism along this path. The de- 
compression following the initial collisional shock should 
be strong enough to push the system inside the unsta- 
ble region of the phase diagram and fragments could be 



formed due to the development of mean-field spinodal 
instabilities 0, 0. However, already in the high den- 
sity phase nucleon correlations are expected to be rather 
large, due to the huge amount of two-body nucleon- 
nucleon collisions. Hence some memory of these high 
density correlations could be kept along the fragment for- 
mation process. 

According to classical molecular dynamics 0, Q or 
lattice-gas calculations [l^, self-bound clusters are ob- 
served, even at equilibrium, in high density systems. This 
has recently suggested an interpretation of the multifrag- 
mentation phenomenon in terms of a sudden explosion 
of the system, where pre-fragments start to appear al- 
ready in the high density phase and they subsequently 
fly apart from each other due to the strong Coulomb re- 
pulsion |16|. However, one should keep in mind that this 
clustering effect revealed in the high density phase could 
be much stronger in classical systems than in nuclear 
matter, which is a Fermi liquid. For instance, mean-field 
calculations of the response of nuclear matter to the pres- 
ence of density fluctuations show that correlations con- 
structed in the high density phase are damped. Indeed, 
according to a mean-field description, it is not energeti- 
cally convenient for a Fermi system at density larger than 
the normal value to develop these high density bumps • 

However, if a large amount of fluctuations is present 
in the high density stage of the reaction, and the system 
expands rather quickly, these correlations could survive 
anyway and play an important role in determining the 
properties of the formed fragments. On the other hand, 
if fluctuations are not so large, these structures could 
be re-organized again by the low density unstable mean- 
field. 

Hence it is clear that the fragmentation mechanism is 
very sensitive to the delicate balance between many-body 
correlations, mean-field effects and the time scale im- 
posed by the reaction dynamics, especially in the Fermi 
energy domain, where one- and two-body effects are 
equally important. Therefore, depending on the way 
one treats the different effects, one could expect a dif- 
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FIG. 1: Contour plots of the density projected on the reaction plane calculated with BGBD for the central reaction 112 Sn + 
112 Sn at 50 MeV/nucleon, at several times (fm/c). The size of each box is 40 fm. 




FIG. 2: The same as in Fig. Qbut calculated with AMD. The size of each box is 80 fm. 



ferent outcome from the available theoretical models de- 
scribing multifragmentation. In this paper we undertake 
this kind of investigation, by comparing the results given 
by two fragmentation models: the stochastic mean-field, 
including momentum dependence (BGBD) |T^ and the 
antisymmetrized molecular dynamics (AMD) models 6]. 
Both have proved to give a good reproduction of some 
aspects of multifragmentation data [l^, 113 ■ We analyze 
central reactions, where we expect to see larger dynam- 
ical (compression-expansion) effects on fragment forma- 
tion. The present study will allow us to get a deeper in- 
sight into the reaction mechanism, in connection with the 
ingredients of the two models. This should be reflected 
into the properties of the obtained primary fragments 
and eventually on measured observables. In this way one 
can also try to identify the experimental quantities that 
are more sensitive to the fragmentation scenario. The 
paper is organized as it follows: In section II we will give 
an outlook of the main ingredients of both models. In 
section III we study the fragmentation path, as given by 
the two models, in the case of a central reaction, 112 Sn 
+ 112 Sn at 50 MeV/nucleon Conclusions are drawn 
in section IV. 



II. INGREDIENTS OF THE MODELS 

Two different kinds of microscopic approaches have 
been proposed and applied to study heavy ion reac- 
tion mechanisms, i.e. to describe the dynamics of nu- 
clear many-body systems. One is the class of molecu- 
lar dynamics models [H EH, EE HI HI while the other 
kind is represented by stochastic mean-field approaches 



num. 

In the mean-field class of descriptions, the dynamical 
state of the nuclear system is characterized by the re- 
duced one-body density in phase space, /(r,p,i), the 
classical analogue of the Wigner transform of the one- 
particle density matrix. At low energies, the time 
evolution of the one-body density is governed by the 
Vlasov equation, which can be regarded as the semi- 
classical approximation to the time-dependent Hartree- 
Fock theory. The residual direct collisions between the 
constituent nucleons are incorporated by means of a 
Pauli-blocked collision integral, leading to the so-called 
Boltzmann-Uehling-Uhlenbeck (BUU) or Boltzmann- 
Nordheim- Vlasov (BNV) approaches [H HI- Tne 
stochastic extension of the transport treatment for the 
one-particle density is obtained by introducing a stochas- 
tic term representing the fluctuating part of the collision 
integral [25j , in close analogy with the Langevin equation 
for a brownian motion. 

Molecular dynamics models usually assume a fixed 
Gaussian shape for the single particle wave functions. 
The many-body state is represented by a simple product 
wave function, with or without antisymmetrization. In 
this way, though single particle wave functions are sup- 
posed to be independent (mean-field approximation), the 
use of localized wave packets induces many-body correla- 
tions both in mean-field propagation and hard two body 
scattering (collision integral), that is treated stochasti- 
cally. From the point of view of stochastic mean-field 
models, the philosophy of molecular dynamics would be 
to introduce a special kind of fluctuation by stochasti- 
cally localizing the single particle wave functions. 

Below we will give the ingredients of the two models, 
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that can be seen as representative of each class, that we 
consider in the present work. 



A. BGBD 

The stochastic mean-field model considered here is 
a semi-classical non-relativistic transport approach, of 
BNV type, see [2^, |3(j , using an isospin and momentum 
dependent effective interaction. The latter is derived via 
an asymmetric extension of the Gale-Bertsch-DasGupta 
(GBD) force EUOH. 

The e nerg y density can be parametrized as follows (see 
also 0,[a|]): 

e = e kin + e(A', A") + s(B', B") + e(C, C") (1) 

where e^in is the usual kinetic energy density and 

e(A',A") = (A 1 + A"f3 2 )^- 
Po 

e(B',B") = (B' + B"{1 2 ) (j\p 

e{C',C") = C'{l NN +lpp) + C"l NP (2) 

The variable (3 = (N—Z) j A defines the isospin content of 
the system, given the number of neutrons (N), protons 
(Z), and the total mass A = N + Z; the quantity po 
is the normal density of nuclear matter. The momentum 
dependence is contained in the T TT i terms, which indicate 
integrals of the form: 

1tt> = J dpdp'f T (r,p)f T ,(r,p')g(p,p') 

where g(p,p') = 1/(1 + (p— p') 2 / A) and f T represents 
the one-body distribution function of neutrons or pro- 
tons. A is a constant, that is taken equal to (1.5 fcp-) 2 , 
being kp the Fermi momentum at normal density. This 
choice of the function <?(p, p') gives a similar behaviour 
with respect to the Gogny effective interaction used in 
the AMD simulations (see next sub-section). We use a 
soft equation of state for symmetric nuclear matter (com- 
pressibility modulus Knm(po) = 215 MeV). In this frame 
we can easily adjust the parameters in order to fix the 
density dependence of effective mass and symmetry en- 
ergy. 

So we describe the time evolution of the system in 
terms of the one-body distribution function / T (r, p, t), 
as ruled by the nuclear mean-field (plus Coulomb inter- 
action for protons) and hard two-body scattering, ac- 
cording to the so-called Boltzmann-Langevin equation 

^ = ^ + {fr,H} = I T [f]+SI T [fl (3) 

where H (r, p, t) is the one-body Hamiltonian, I T [f] is the 
average two-body collision integral and SI T [/] represents 



the stochastic source term 0, H^, |2(| . The test particle 
method is used to solve numerically Eq. The free 
energy- and angle-dependent nucleon-nucleon cross sec- 
tion is used in the collision integral. Fluctuations are 
introduced within this mean- field treatment, according 
to the approach presented in [2!||35j], i.e. by agitating 
the spatial density profile. Once local thermal equilib- 
rium is reached, the density fluctuation amplitude, a p 
is evaluated by projecting on the coordinate space the 
kinetic equilibrium value of a Fermi gas. Then, in the 
cell of r space being considered, the density fluctuation 
bp is selected randomly according to the gaussian distri- 
bution exp(— 8p 2 /2a 2 ). This determines the variation of 
the number of particles contained in the cell. A few left- 
over particles are randomly distributed again to ensure 
the conservation of mass. Momenta of all particles are 
finally slightly shifted to ensure momentum and energy 
conservation. Hence, while the dynamical evolution of 
the system is still described in terms of the one-body dis- 
tribution function, this function experiences a stochastic 
evolution, in response to the action of a fluctuation term 
essentially determined by the degree of thermal agitation 
present in the system. 

According to this stochastic mean-field theory, the 
fragmentation process is dominated by the growth of vol- 
ume (spinodal) and surface instabilities encountered dur- 
ing the expansion phase of the considered excited systems 
J36j. Therefore density fluctuations provide the seeds of 
the formation of fragments, whose characteristics are re- 
lated to the properties of the most unstable collective 
modes of the mean-field. In finite nuclei, several multi- 
poles are excited with close probabilities. Hence a large 
variety of fragment configurations may be obtained, due 
to the beating of the several unstable modes 0, Hi| . 
This description of the fragmentation path can explain 
several features, concerning also rather exclusive observ- 
ables 0], observed in experimental data at around 30 
MeV/nucleon 



B. AMD 

For the AMD [T3, |U |37| calculations pre sented here, 
we use the same framework as in Ref. 18] which can 
reproduce the fragment charge distribution of the central 
Xc + Sn collisions at 50 MeV/nucleon. 

In AMD, we employ the Slater determinant of Gaus- 
sian wave packets 



,r A \${Z)) oc det 



exp 



{-i/fa-Z^/^ 2 }*^*) 



(4) 

where Xai are the spin-isospin states with on = p "f, 
p I, n |, or n j. Thus the many-body state \®{Z)) 
is parametrized by a set of complex variables Z = 
{Zj}j = i a, where A is the number of nucleons in the 
system. The width parameter v = (2.5 fm)~ 2 is treated 
as a constant parameter common to all the wave packets. 
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If we ignore the antisymmetrization effect, the real part 
of Zj corresponds to the position centroid and the imag- 
inary part corresponds to the momentum centroid. This 
choice of wave functions is suitable to describe fragmen- 
tation channels where each single particle wave function 
should be localized within a fragment. 

The dynamics of fragmentation is a highly complicated 
quantum many-body problem in which huge number of 
fragmentation channels will appear in the course of the 
evolution. An AMD wave function [Eq. Q] is intended to 
describe one of the channels rather than the total many- 
body state, and the emergence of channels is represented 
approximately by some stochastic terms in the equation 
of motion. 

The stochastic equation of motion for the wave packet 
centroids Z may be symbolically written as 

^ Z, = {Zi , H} PB + (NN coll) +AZ i {t)+fi(Z i ,H'). (5) 
at 

The first term {Zi,7i}pB is the deterministic term de- 
rived from the time-dependent variational principle with 
an assumed effective interaction such as the Gogny inter- 
action j^]. The second term represents the effect of the 
stochastic two-nucleon collision process. The collisions 
are performed with the "physical nucleon coordinates" 
that take account of the antisymmetrization effects, and 
then the Pauli blocking in the final state is automatically 
introduced 22]. The third term AZj(t) is a stochastic 
fluctuation term that has been introduced in order to re- 
spect the change of the width and shape of the single 
particle distribution [(i fl8L l39l | . In other words, the com- 
bination {Zi,7Y}pB + AZj(f) approximately reproduces 
the prediction by mean field theories for the ensemble- 
averaged single-particle distribution, while each nucleon 
is localized in phase space for each channel. The term 
AZj(t) is calculated practically by solving the Vlasov 
equation with the same effective interaction as for the 
term {Zi,7i}pB- In the present version of AMD [l^, the 
property of the fluctuation AZj(i) is chosen in such a 
way that the coherent single particle motion in the mean 
field is respected for some time interval until the nucleon 
collides another nucleon. The last term /x(Zj,W) is a 
dissipation term related to the fluctuation term AZi(t). 
The dissipation term is necessary in order to restore the 
conservation of energy that is violated by the fluctuation 
term. The coefficient fi is given by the condition of energy 
conservation. However, the form of this term is somehow 
arbitrary. We shift the variables Z to the direction of the 
gradient of the energy expectation value TL under the con- 
straints of conserved quantities (the center-of-mass vari- 
ables and the total angular momentum) and global one- 
body quantities (monopole and quadrupole moments in 
coordinate and momentum spaces). A complete formu- 
lation of AMD can be found in Refs. ^Il37|. 

In the present work, we use the Gogny effective inter- 
action [33] which corresponds to a soft equation of state 
of symmetric nuclear matter with the incompressibility 
^nm(po) = 228 MeV. The mean field for this force has 



a momentum dependence which is similar to that in the 
BGBD calculation. As the two-nucleon collision cross 
sections (a pp = a nn and cr pn ), we use the energy- and 
angle-dependent values in the free space with the maxi- 
mum cut-off of 150 mb. In order to avoid low energy spu- 
rious collisions, due to Pauli blocking violation (caused 
by the finite number of test particles), in the BGBD cal- 
culations a lower cut-off, 50 mb, has been considered. We 
have checked that this leads to approximately the same 
degree of stopping in the two models. 



C. Some remarks 

We would like to stress here the main differences be- 
tween the two models described above. 

In the AMD version used here a special proce- 
dure is adopted, that ensures a coherent single-particle 
motion in the mean-field (including diffusion and shrink- 
ing effects of the nucleon wave packet in phase space) 
until the considered nucleon collides another one. It can 
be demonstrated that this procedure reproduces exactly 
the coherent time evolution of the Wigner function / in 
the case of a harmonic oscillator potential (or for free 
nucleons), though this correspondence is not exact in the 
general case. Hence one-body effects should be similarly 
treated in the two models. From this point of view, the 
present version of AMD is rather different from earlier 
molecular dynamics formulations [20L l22j , where the use 
of localized wave packets in the full dynamics implies 
that one-body effects are not as precisely described as in 
mean-field models. 

The most relevant difference between AMD and BGBD 
is related to the method followed to implement stochastic 
two-body scattering. In fact, in BGBD fluctuations are 
introduced by agitating the one body density function, to 
account for the stochastic part of the collision integral, 
according to the Boltzmann-Langevin theory 0, • In 
some sense, this would correspond to a description of the 
system in terms of unrestricted fluctuating single particle 
wave functions. On the other side, in AMD fluctuations 
are introduced by stochastically localizing the single par- 
ticle wave functions in phase space when a two-nucleon 
collision takes place. 

It is difficult to discuss the general validity of the vari- 
ous approximations adopted to solve the quantum many- 
body dynamics because it may depend on the particular 
reaction mechanism and energy range under study. Here 
we will focus on the description of the multifragmenta- 
tion mechanism at 50 MeV/nucleon. While it is well 
known that the results of early molecular dynamics mod- 
els (23, H3| and standard mean-field approaches [2^ are 
rather different, the improved AMD and the stochastic 
BGBD can be considered as closer approaches and it is 
interesting to investigate how their respective predictions 
compare to each other in the case of fragmentation reac- 
tions. 
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III. ANALYSIS OF THE RESULTS 

Our aim is to investigate the fragmentation path in 
violent collisions at intermediate energy, as predicted by 
the BGBD and the AMD models. As discussed above, in 
the BGBD model one essentially follows the evolution of 
the one-body density and the fragmentation mechanism 
is mainly based on the amplification of its fluctuations. In 
AMD nucleon wave packets are propagated, from which, 
however, it is possible to reconstruct the one-body den- 
sity and related observables. 

Hence our study of the reaction path will be per- 
formed by looking at quantities connected to the one- 
body density and to its fluctuations. We will investi- 
gate central collisions of the system 112 Sn + 112 Sn at 
50 MeV/nucleon. An ensemble of 200 trajectories with 
b = 0.5 fm has been collected with BGBD, while in AMD 
20 events with < b < 1 fm are considered. 



A. One-body observables 

To give a qualitative representation of the time evolu- 
tion of the system, density contour plots in the reaction 
plane, as obtained in the two models, are shown in Figs, 
n and |21 at several time steps. As one can see from the 
figures, both models predict that the system is initially 
compressed. Then expansion follows and several inter- 
mediate mass fragments (IMF) appear. 

In BGBD, according to the value of the spatial den- 
sity p(r), one can identify a "gas" phase (p < po/3), 
associated with particles that leave rapidly the system 
(pre-equilibrium emission) and/or are evaporated, and a 
"liquid" phase, where fragments belong to. po denotes 
the normal density value. As a matter of fact, we ob- 
serve that the "liquid" essentially corresponds to parti- 
cles with charge greater than 2. In AMD, the fragments 
with A < 4 are regarded as in the gas part. The time evo- 
lution of the number of nucleons that are in the gas phase 
is represented in Fig. [21 It is possible to observe that the 
number of particles that escape at early times from the 
interacting nuclear matter, due to two-body scattering 
(pre-equilibrium effects), is different in the two models. 
In the BGBD case (full line) the emission rate is larger 
with respect to the AMD calculations. Moreover, there 
is a clear change of slope at around t w 100 fm/c where, 
as we will see in the following (see also Fig.l), the nuclear 
system reaches low density values, fragments start to be 
formed and nucleons are emitted due essentially to evapo- 
rative processes. This is not so evident in the AMD case. 
This result seems to indicate already different emission 
mechanisms between the two models. In BGBD the pro- 
duction mechanism of pre-equilibrium light particles and 
the following evolution of the system (fragmentation) are 
somehow decoupled, while in the AMD case there are no 
easily identifiable emitting sources and all fragments and 
light particles seem to be emitted on about equal footing. 
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FIG. 3: Time evolution of the number of nucleons in the "gas" 
phase, as obtained in the central reaction 112 Sn + 112 Sn at 50 
MeV/nucleon. Full line refers to BGBD calculations, while 
dashes represent AMD calculations. 



To closely follow the time evolution of the system and 
to better characterize its fragmentation path, we have 
studied the behaviour of the following observables: the 
radial density profile of the system and the radial collec- 
tive momentum as functions of time. The radial density 
at a given distance r is obtained by averaging the lo- 
cal density p(r) over the surface of a sphere of radius r. 
The radial collective momentum is the projection of the 
collective momentum at the position r along the radial 
direction, averaged over the surface of the sphere of ra- 
dius r. These quantities are further averaged over the 
event ensemble. 

In the BGBD case, the behaviour of the radial den- 
sity profile, presented in Fig. 21 indicates that, after an 
initial compression (t = 40 fm/c), the system expands 
and finally it gets rather dilute, due to the occurrence 
of a monopole expansion, generated by the compression. 
As one can see from Fig. ^ while the system expands, it 
breaks up into many pieces. The matter appears mostly 
concentrated within a given window of the radial distance 
(see for instance the line at t = 120 fm/c), indicating the 
formation of a bubble-like configuration, where fragments 
are located. Indeed the central region of the system is 
rapidly depleted. 

It is rather interesting to look at the profile of the col- 
lective momentum (divided by the nucleon mass), v(r), 
see Fig. The time evolution of this quantity is largely 
influenced by the occurrence of monopole compression 
and expansion. Indeed, after the initial compression, the 
restoring force generated by the mean-field pushes the 
system back to the normal density value and the collec- 
tive momentum increases (compare the evolution from 
t = 40 to 60 fm/c). It can be noticed that, at t = 60 
fm/c, in the region where the "liquid" is located, i.e. for 
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FIG. 4: Density profiles, at several times (in fm/c), as ob- 
tained in the BGBD case. 



FIG. 6: Density profiles, at several times (in fm/c), as ob- 
tained in the AMD case. 
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FIG. 5: Collective momentum profiles, in units of the nucleon 
free mass, at several times (in fm/c), as obtained in the BGBD 
case. 



FIG. 7: Collective momentum profiles, in units of the nucleon 
free mass, at several times (in fm/c), as obtained in the AMD 
case. 



r < 8 - 9 fm (cf. Fig.QJ), the collective momentum profile 
is almost self-similar, i.e. the radial momentum is propor- 
tional to the radial distance. Then the system goes into 
the expansion phase and the collective momentum de- 
creases again. This deceleration indicates that the system 
is still rather homogeneous while it expands, reaching the 
low density (unstable) region. Then fragmentation can 
be associated with the occurrence of spinodal decomposi- 
tion ■ More specifically, the system slows down due to 
the presence of a counter-streaming flow that develops 
from the surface towards the interior, trying to recom- 
pact the system. This acts against the initial expansion, 
so the collective momentum profile is modified and, in 
the region where the matter is located (cf. Fig. ^} , the 
system is slowed down. So actually when fragments start 
to appear (t rs 100 fm/c) their collective momentum is 
not so large. However, this is rather large in the central 
part of the system (r < 3 fm at t = 60 and 100 fm/c), 
which is related to the rapid depletion of the density in 
this region. 

One may also notice that, at large radial distance, 
where mostly pre-equilibrium particles are located, the 
radial momentum exhibits a different trend and is rapidly 
increasing with r. 



In AMD calculations, the density profile in Fig.[fJ]shows 
the time evolution of the compression and expansion 
which is qualitatively similar to the BGBD case. How- 
ever, we notice that AMD shows broader density distribu- 
tion than BGBD as the system expands. The bubble-like 
configuration is not as prominent as in BGBD. The anal- 
ysis of the radial momentum profile (Fig. 0) also reveals 
that the expansion dynamics in AMD is different from 
the BGBD case. When the system goes from the com- 
pressed state towards the normal density value (i « 70 
fm/c), the collective momentum increases less than in 
the BGBD case. This suggests that the energy that was 
stored in compression goes into larger (with respect to 
BGBD) thermal kinetic fluctuations, which is consistent 
with the broader density profile observed in Fig. El It 
should be noticed that the broad density profile in the 
late stage does not mean homogeneous dilute matter but 
corresponds to the situation in which fragments are dis- 
tributed widely in the space. By the comparison with 
the BGBD case, it seems that the system ceases to be- 
have as homogeneous matter already at around t = 70 
fm/c, before entering the low density (spinodal) region. 
Then, contrarily to what happens in the BGBD case, 
where the system is slowed down during the expansion 
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FIG. 8: Time evolution of the density variance, as obtained FIG - 9: Collective momentum variance, as obtained in BGBD 
in BGBD (full line) and AMD (dashed line) calculations. ( ful1 line ) and AMD (dashed line) calculations, as a function 

of the distance on the z axis. 



phase, after t = 70 fm/c the expansion collective momen- 
tum keeps almost unchanged in the AMD calculations. 
The mean field restoring force, that would recompact the 
system, appears less effective in the event-averaged one- 
body dynamics. Moreover, after this time, the radial 
dependence of the collective momentum is almost self- 
similar for the 'liquid' part. At all the time, we see a 
kink in the v(r) curve. The exterior and rapid compo- 
nent corresponds to the pre-equilibrium particles and the 
inner and slower component corresponds to the 'liquid' 
composed mainly of fragments. However the kink is al- 
ways more pronounced in the BGBD case, pointing to- 
wards a larger difference of velocity between fragments 
and pre-equilibrium particles. 



B. Fluctuations 

The analysis of the density variance is an important 
tool to identify the moment when fragmentation sets in 
and to get some information about the related mecha- 
nism. During the fragmentation process, the density vari- 
ance, evaluated at a given position r, grows in time, then 
it saturates (when fragments do not interact anymore, 
apart from the Coulomb repulsion) and eventually it de- 
creases due to the fact that fragments fly apart from each 
other. For instance, in the spinodal decomposition sce- 
nario, the density fluctuations grow exponentially, with a 
characteristic amplification time r (typically of the order 
of 30-50 fm/c) related to specific properties of the nuclear 
interaction such as its range Q- 

Density variances have been calculated starting from 
the value of the density along the x-, y- and z-axes. For 
instance, the variance of the density along the z-axis is 
defined by S%(z) = ((p(0,0,z) - (p(0, 0, z))) 2 ), where 
p(x, y, z) is the density for each event and the brackets 
( ) stand for the average value over the event ensemble. 

In Fig. [S] we show the time evolution of S z , which is a 



representative value of S z {z), defined by 

_ JS z (z)(p(0,0,z))dz 

j{p(0,0,z))dz ■ {0) 

Similar behaviour is observed for the variance in the other 
directions, that can be defined analogously. One can see 
that density fluctuations have close values in the two ap- 
proaches in the high density phase (t « 40 fm/c). How- 
ever, in the BGBD case fluctuations are damped while 
the system expands, relaxing towards the lower equi- 
librium value expected for nuclear matter at lower den- 
sity and temperature. Only when low-density values are 
reached (t w 70 fm/c) and the mean-field is unstable, 
one starts to see the rapid increase of density fluctua- 
tions, corresponding to fragment formation. This is due 
to the fact that density bumps are amplified, leading to 
fragments at normal density. When fragments fly apart 
from each other at large time t, the value of S z gradually 
decreases due to a trivial effect of the definition of S z 
[Eq. ©]. 

On the other hand, in the AMD case the density fluc- 
tuations continue to grow as the system expands from 
t « 40 to 70 fm/c, suggesting that prefragments start 
to develop gradually at this stage even though the den- 
sity is not very low. This is probably due to the nucleon 
wave packet localization when two-body scattering oc- 
curs. As we have already mentioned, the relatively small 
magnitude of the collective momentum observed in Fig.0 
around t sa 70 fm/c may be associated with the large mo- 
mentum fluctuation. If the velocities of prefragments are 
randomly distributed, the collective momentum will fluc- 
tuate from one event to another. This is confirmed by the 
analysis of collective momentum fluctuations, displayed 
in Fig. EI In i&ct, when the system is expanding (t = 60 
- 70 fm/c), larger fluctuations are seen in the AMD case. 
Furthermore, the mean-field restoring force for the global 
expansion dynamics is expected to be less effective in this 
situation and hence deceleration of expansion is weak in 
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the AMD case as observed in the comparison of Figs. [5] 
andQ 

IV. DISCUSSION 

According to the results presented above, emission 
mechanisms appear rather different in the two mod- 
els. For light particles, a more abundant pre-equilibrium 
emission (Z < 3) is observed in the BGBD case. Con- 
cerning the IMF emission mechanism, the involved time 
scales and the relative importance of one- and many-body 
effects appear different in the two models, due to the dif- 
ferences in the implementation of mean-field propagation 
and two-body scattering. In the BGBD case the fragmen- 
tation process follows the spinodal decomposition sce- 
nario. Indeed the system enters the low density region as 
a nearly homogeneous source, then density fluctuations 
are amplified. The formation of bubble-like structures is 
favoured, as a manifestation of monopole instabilities. In 
the AMD model, many-body correlations have a stronger 
impact on the fragmentation dynamics, while mean-field 
effects appear weaker. Fragment formation seems not 
to be directly linked with the spinodal decomposition 
scenario, but rather with the appearance of some pre- 
fragments at earlier times. As a consequence, fragments 
are formed on shorter time-scales in AMD. Finally, we 
observe that pre-equilibrium and IMF emissions happen 
on different time scales and are related to different mech- 
anisms in the BGBD case, while in AMD the distinction 
of the two time scales is not as clear as in BGBD. 

As shown in Fig. [8] spatial density fluctuations have 
similar values, in the high density phase, in the two mod- 
els. However, when the system starts to expand, fluctu- 
ations are quenched in BGBD (until mean-field insta- 
bilities are encountered) while they increase in AMD. 
The same is true for the collective momentum fluctua- 



tions (see Fig. |5J. These differences are naturally un- 
derstood as an effect of the nucleon localization in the 
AMD case, that induces larger many-body correlations. 
One can conclude that, in the energy range considered 
here (50 MeV/nucleon), the fragmentation path is really 
sensitive to the interplay between one- and many-body 
effects. Changing the relative weight of these effects leads 
to a rather different outcome. 

As a consequence of the emission mechanisms outlined 
above, one expects to see different primary fragment con- 
figurations in the two models, with larger fluctuations in 
the AMD case, though results could be similar for inclu- 
sive observables, such as charge distributions, especially 
after secondary decay effects have been considered. Con- 
sidering kinematic properties and more exclusive observ- 
ables, such as IMF-IMF correlations, that keep a better 
track of the freeze-out configuration and primary frag- 
ment partitions, one should be able to disentangle be- 
tween the predictions of the two models. A compari- 
son with available experimental data would allow to shed 
some light on the mechanisms and most relevant effects 
that govern the fragmentation process. This analysis will 
be the subject of a forthcoming paper. 

Also, it would be interesting to extend this study 
to semi-peripheral collisions, to investigate transport 
properties and dissipation-equilibration mechanims. In 
charge asymmetric reactions, isospin equilibration could 
be a good tracer of the reaction dynamics. 
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